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Abstract 

The paper studies the effects of a time-delayed feedback control on the appearance and 
development of spatio-temporal patterns in a reaction-diffusion system. Different types of control 
schemes are investigated, including single-species, diagonal and mixed control. This approach 
helps to unveil different dynamical regimes, which arise from chaotic state or from tavelling 
waves. In the case of spatio-temporal chaos, the control can either stabilize uniform steady 
states or lead to bistability between a trivial steady state and a propagating travelling wave. 
Furthermore, when the basic state is a stable travelling pulse, the control is able to advance 
stationary Turing patterns, or yield the above-mentioned bistability regime. In each case, the 
stability boundary is found in the parameter space of the control strength and the time delay, 
and numerical simulations suggest that diagonal control fails to control the spatio-temporal 
chaos. 

Various real-life systems exhibit complex dynamical regimes, including spatial pat- 
terns and spatio-temporal chaos. An important issue is the possibility of attaining 
a desired state by means of some external influence on the system. Time-delayed 
feedback uses the difference between current state of the system and its state some 
time ago to provide an efficient tool to control system's dynamics. In particular, this 
approach allows one to effectively switch between different spatio-temporal regimes 
in the system. In this paper, we show how the time-delayed feedback control can be 
used to control spatio-temporal chaos, as well as to provide a transition from travelling 
waves to stationary periodic patterns. 
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1 Introduction 



Spatially extended systems arise in modelling of various physical, biological, chemical and engineer- 
ing phenomena, and an ultimate understanding of such systems gives one the ability to control them 
in order to achieve desired behaviour or spatial patterns. Since many systems in nature exhibit 
intrinsic chaotic behaviour, control of chaos in spatially extended systems has recently become an 
active and important area of research. 

Control of a dynamical system in a chaotic state usually requires application of some external 
perturbation in order to achieve a desired type of behaviour, such as a steady state, regular or 
quasi-periodic oscillations. A number of different methods and techniques have been suggested 
in order to control chaos pp. For systems without spatial extension, standard approaches include 
the Ott-Grebogi-Yorke (OGY) scheme [2] and time-delayed feedback control (TDFC) [3]. TDFC 
schemes are straightforward to implement and can also be adapted to control spatially extended 
systems. These schemes are based on the use of the difference between the system variables at the 
current moment of time and their values at some time in the past. For low-dimensional systems, 
TDFC has been successfully used to stabilize unstable steady states [HI IH E], as well as unstable 
periodic orbits [7] in a number of physical and biological systems [I]. 

More recently, TDFC has been used to control spatio-temporal chaos and pattern formation in 
chemical reactions. Close to the onset of uniform oscillations, it is possible to reduce any reaction- 
diffusion system to an amplitude equation, called the complex Ginzburg-Landau equation (CGLE) 
[8]. The CGLE helps to understand the underlying dynamics of the reaction-diffusion systems 
close to the bifurcation point, however, it does not provide a lot of information about the dynamics 
far away from the onset of oscillations. Beta and Mikhailov have shown analytically that global 
TDFC cannot stabilize uniform oscillations when the system in the spatio-temporally chaotic regime 
[10] ■ At the same time, local time-delay feedback control which also includes spatial feedback 
terms has been shown to be effective in stabilizing travelling waves in CGLE in one and two spatial 
dimensions . 

In globally coupled reaction-diffusion systems, which arise in modelling semiconductor nanos- 
tructures, global or local TDFC and its various modifications have proved successful in stabilizing 
unstable spatio-temporal oscillations [131 El US]- Unstable rigid rotation of spiral waves has also 
been stabilized by TDFC |16j . Pulse propagation in the FitzHugh-Nagumo reaction-diffusion model 
has recently been shown to be controllable by local TDFC or nonlocal instantaneous feedback [17J. 

From a practical perspective, there are several issues that have to be taken into account when 
considering implementation of specific TDFC schemes. One such aspect concerns the spatial range 
over which the control has to be applied. Ahlborn and Parlitz [18J have recently shown how one can 
use multiple time-delayed control applied at a number of control points to stabilize plane waves in 
the CGLE, while Bleich and Socolar [19] have illustrated stabilization of periodic orbits in chaotic 
CGLE using local control. Another problem is that most examples of TDFC for spatially-extended 
systems such as, for instance, one-dimensional chaotic CGLE [19] and spatio-temporally chaotic 
semiconductor laser arrays |2Uj . use diagonal control with symmetric gain, which may be quite 
restrictive. In the case of optical systems such as laser devices it is possible to have a non-diagonal 
coupling, realized experimentally using an optical phase as an additional control parameter [HE]. 
Hence, it is important to consider the influence of various TDFC schemes on the dynamics of the 
system and analyze different types of spatio-temporal dynamical regimes that can be achieved. 

In this paper we perform a systematic study of the effects of local TDFC on spatio-temporal 
dynamics in a two-component reaction-diffusion system. As a testbed for this analysis we use the 
paradigmatic Gray-Scott model derived in the context of chemical reactions. It has already been 
shown using extensive numerical simulations that the Gray-Scott model supports a wide range 
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of spatio-temporal dynamics, including stationary inhomogeneous patterns, travelling fronts and 
pulses, as well as spatio-temporal chaos. A particular feature of this system, which makes it different 
from other reaction-diffusion models, is the occurrence of pulse-splitting, where a travelling pulse 
leaves in its wake other pulses propagating in different directions. 

The outline of this paper is as follows. In the next section we introduce the Gray-Scott system. 
Section 3 contains an analysis of the effects of TDFC on spatio-temporal chaos. In Section 4 various 
TDFC schemes are applied to propagating pulses. The paper concludes with a summary of our 
findings, together with a consideration of their implications. 



2 Gray-Scott system 

In order to analyse the effects of TDFC on the spatio-temporal dynamics in reaction-diffusion 
systems, we consider the Gray-Scott model, which describes a cubic auto-catalytic chemical reaction 
of the type [21] 

U + 2V^3V, 

V -> P, { ' 

where U is continuously supplied into an open flow reactor, and the product P is removed. The 
first equation represents an autocatalytic process, in which two molecules of species V produce 
3V through interaction with a molecule of species U. If the system were closed, then by virtue 
of irreversibility of the two reactions, all reactants would eventually turn into a product and be 
removed. However, when the reactor is continuously fed with a uniform flow of species U, it is 
possible to maintain far-from-equilibrium conditions, in which a wide range of possible dynamics 
can be observed. In experiments, a gel is used as a medium to prevent the occurrence of convective 
currents [22]. The kinetic equations for the system ([!]) can be written after some rescaling as 

— = -uv 2 + o(l -u) + D u V 2 u = f(u, v) + D u V 2 u, 

(2) 

— = uv 2 - (a + b)v + D v V 2 v = g(u, v) + D V V v, 

where u and v are the concentrations of the species U and V, respectively, a is the inflow rate, a + b 
is the removal rate of V from the reaction, and D u and D v are the diffusion coefficients of the two 
species. 

The system ^ can possess up to three homogeneous steady states. There is a trivial steady 
state 

E o = (u o ,v o ) = (l,0), (3) 
which exists for all parameter values and corresponds to a steady flow without reaction. Provided 

d = 1 -A{a + b) 2 /a > 0, (4) 

the system has two non-trivial steady states 

E 1 = (u 1 ,v 1 )= Q(l-V3),i^(l + V^)), (5) 

and 

E2 = (u 2 ,v 2 )= Q(l + V2), —^(1-^/3)). (6) 
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Figure 1: (a) Bifurcation diagram of the Gray-Scott system ([2]) without diffusion. Solid lines 
correspond to stable steady states, dashed to the unstable steady state E2, dotted to the Hopf 
unstable steady state E\, and dash-dotted indicates a saddle- node (SN) line. The cross (X) indicates 
the location of the Takens-Bogdanov point, (b) Phase diagram in (a, b) space. Solid line is the 
boundary of saddle-node bifurcation, dotted line is the boundary of Hopf bifurcation, and dash- 
dotted line is Turing instability boundary. TB is the Takens-Bogdanov point. 



Linear stability analysis indicates that the trivial steady state Eq is always linearly stable (even 
with respect to spatially inhomogeneous perturbations), while the steady state E2 is unstable with 
respect to homogeneous perturbations for any parameter values for which it exists. Figure [T] shows a 
bifurcation diagram for the spatially homogeneous Gray-Scott model, as well as a phase diagram in 
(a, b) space. One of the features is that when the inflow a is large enough (a > 1/4), the production 
of V is suppressed, and as a result the only fixed point of the system is the trivial steady state Eq, 
as shown in Fig. [ija). 

In Fig. [T|(b) the saddle-node bifurcation curve is given by 



a 



SN 



1 



[l-86± y/1- 166]. (7) 



To the right of a only the trivial steady state Eq exists. As a SN is crossed, two more states E\ 
and E2 appear. While the state E2 is always unstable, the fixed point E\ may lose its stability via 
a Hopf bifurcation. The boundary of this bifurcation is given by 



Vb-2b- Jb(l-4Vb) 



(8) 



This boundary starts at (a, b) = (0, 0) and follows the lower branch of the saddle-node curve until 
it coincides with it at a Takens-Bogdanov point (a, b) = (1/16, 1/16). The steady state E\ is stable 
above the Hopf boundary and unstable below it. It is noteworthy that the stability of the limit 
cycle that appears due to the Hopf bifurcation varies along the Hopf curve: it is stable for b < 0.035 
and unstable for b > 0.035. 
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Another way in which a state E\ may lose its stability is via a Turing bifurcation, which may 
give rise to stationary inhomogeneous patterns. Linear stability analysis in this case shows that 
the Turing bifurcation occurs when the ratio of the two diffusion coefficients a = D u /D v satisfies 
the condition [23j 

[a(a + b) - (vf + a)] 2 > 4a(a + b)(vf - a), (9) 

where v\ is the steady-state value of the species V given in equation ([5]). The higher a, the wider is 
the region in (a, b) space where Turing bifurcation can be observed. The steady state E\ is stable 
with respect to spatially inhomogeneous perturbations above the Turing boundary and unstable 
below it. For sufficiently high values of cr, the boundary of Turing bifurcation will move above 
the boundary of Hopf bifurcation. For a = 2, which will be used in numerical simulations later 
in this paper, the region of Turing instability is almost entirely contained within the parameter 
region of Hopf instability, as shown in Fig. Qb). Interaction of these two instabilities gives rise 
to spatio-temporal chaos and mixed modes having time-dependent spatial structures |23| 124] . In 
Figure Qb), the region of spatio-temporal chaos is narrow and bounded by the Turing boundary 
from the top and the saddle-node curve from the bottom. It is noteworthy that very close to the 
Takens-Bogdanov point one has a Turing instability where stationary inhomogeneous patterns are 
observed. 

It has already been mentioned that the Gray-Scott system exhibits a wide range of dynamical 
behaviours. In the parameter regime away from Turing and Hopf areas, the system always ends up 
in either the homogeneous trivial steady state Eq or state E\. The simplest inhomogeneous state is 
represented by a stationary non-uniform pattern arising from a Turing bifurcation. Besides these 
temporally stationary structures, the Gray-Scott system is able to support standing or travelling 
waves and fronts (the latter refer to heteroclinic connections between Eq and E\). Travelling waves 
occur in the parameter region outside the saddle-node boundary. At a higher level of complexity 
are the mixed modes, describing time-dependent spatial structures, the so-called self-replicating 
patterns and spatio-temporal chaos |24j . The self-replicating patterns represent heteroclinic con- 
nections between the trivial homogeneous state and a spatially periodic pattern. 



3 Control of spatio-temporal chaos 

In this section we concentrate on the parameter regime where the system exhibits spatio-temporal 
chaos and consider several schemes for the control of this state based on TDFC. This approach has 
already been used to control spatio-temporal chaos in globally coupled reaction-diffusion systems 
\l'S\ [T5] , as well as in reaction-diffusion systems with global control [9] ES] • 
Formally, the locally controlled Gray-Scott model can be written as 

d f u\ _ f f(u,v)\ f D u \ 2 f u\ f u(t-r)-u(t)\ 

dt{v )-{ g (u,v) ) + {o D v ) U ) + KA {v(t-r)-v(t) J' (10) 

where K is the control strength which can be either positive or negative, r > is the time delay, 
and A is a two-by-two matrix which describes the particular coupling of the control term. First we 
consider the case when the control of each variable depends only on the history of that variable. 
For single-species control, the matrix A takes the form 

A u=(i !l or A v =( ° ° ), (11) 



x J ' " \ 1 / 

for control only in u (activator control) and only in v (inhibitor control), respectively. 
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Figure 2: (Color online) Stability boundaries of the homogeneous steady state E± from the char- 



acteristic equation (12). In both cases, the steady state E\ is stable above the boundary for q = 
described by (14). (a) Activator control, (b) Inhibitor control. Parameter values are a = 0.028, 
b = 0.053, D u = 2- 10" 5 , D v = 10~ 5 . 



In the case of activator control, we look for solutions of the system ( 10 ) linearized near E\ 
(ui, v\) in the form 



C 2 



D iqx+\t 



where u = u\ + u, v = v\ + v, C\ i are constants, q is the wavenumber of the perturbation, and A 



is the corresponding growth rate. Substituting this into the system (10) with A u from (11) yields 
the following characteristic equation 



D u q 2 + K[e 



1]-A -2(o + 6) 
a + b — D v q 2 - 



X 



0. 



(12) 



In the absence of control (K = 0), it follows from this equation that the uniform mode with q = 
is the most unstable, i.e. it has an eigenvalue with the largest positive real part. Since the diffusion 
plays a stabilizing role, all other modes with q > will be more stable, as is confirmed by numerical 
solution of equation (12) shown in Fig. [2] For this reason, it suffices to control the uniform mode 
q = to achieve complete stabilization of the steady state E\ for any spatially inhomogeneous 
perturbation. When q = 0, the equation (fl2|) can be written as 



A 2 + A 



b-K e 



-At 



+ (a + b) 



a + K(e 



-At 



0. 



(13) 



In order to find the control strength K that can stabilize the steady state E\ for each value of the 
time delay t, we compute the stability boundary in the (K, r) plane by setting A = iuj in equation 



(13). Separating real and imaginary parts gives the following system 



-ui 



+ (a + &)(«? 



K) — Ku sin lot + {a + b)K cos cot = 0, 



cj(— b + v I + K) — ujK cosujt — (a + b)K sinwr = 0. 
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Figure 3: (Color online) Control domain of spatio-temporal chaos, (a) Activator control, (b) 
Inhibitor control, (c) Phase diagram for non-diagonal control (15). Colour denotes a specific final 



state: spatio-temporal chaos (black), stable mixed (Turing-Hopf) state (pink), coarsening (green), 
bi-stability between travelling waves and a trivial steady state E$ (red) and a uniform non-trivial 
steady state E\ (yellow). Parameter values are a = 0.028, b = 0.053, D u = 2 • 10 -5 , D v = 10~ 5 . 



This system can be solved to yield the pair (K, r) as parameterized by oj: 



K 



1 [(q + b) 2 + oj 2 ] [(a - v 2 ) 2 + a; 2 ] - Auj 2 bv 2 

2 (a- Vl ) 2 (a + b) 2 + u 2 {a + v 2 ) ' 



(14) 



1 



Arct 



an 



2oj 4(a + b)vfoj 
~K ~ K [(a + b) 2 + uj 2 ] 



db mr 



n 



0,1,2,.., 



where Arctan corresponds to the principal value of arctan. We are only interested in the branch 
n = 0, as for other values of n the corresponding control term in the system ( 10 ) will be very large 



compared to other terms, making it unfeasible for practical purposes. In Fig. [2^ a) we show this 
stability boundary, which corresponds to q = 0, together with stability boundaries when spatially 
inhomogeneous perturbations are taken into consideration. For any q > 0, these boundaries lie 
below the stability boundary for the uniform mode, as we have explained earlier. Figure [2jb) shows 
the results of computation of the stability boundary in the (K, r) plane for inhibitor control. One 
can note the approximate symmetry between the stability boundaries for activator and inhibitor 
control up to a reverse of the sign of the control strength. This is due to the fact that the inhibitor 
control would result in the characteristic equation similar to (13), but with the term (a+b)K[e~ XT — 
1] being replaced by — (v \ + a)K [e~ Ar — 1] , where the difference in the absolute values of (a + b) and 
{v\ + a) is less than 10%. Finally, we note that it is only possible to stabilize the steady state E\ 
for negative values of the control strength K in the case of activator control, and only for positive 
K in the case of inhibitor control. 



In order to understand the full dynamics of the system ( 10 ) with time-delayed feedback control 



we solve this system numerically and record the ultimate stage of time evolution. System parameters 
are fixed and taken to be D u = 2 • 10~ 5 and D v = 10~ 5 , while kinetic parameters are a = 0.028 
and b = 0.053. In each case we solve system (10) on an interval of length L = 2.5 using an explicit 
Euler scheme with a spatial resolution of 256 mesh points, a time-step of At = 0.05 and periodic 
boundary conditions. Initial conditions are chosen to be a random localized perturbation of Eq 
in the middle of the domain which then evolves into a spatio-temporally chaotic state. After the 
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Figure 4: (Color online) Space-time plot in the case of activator control, (a) Spatio-temporal chaos, 
(b) Non-trivial steady state E\. The colour code corresponds to the values of u(x,t). Parameter 
values are a = 0.028, b = 0.053, D u = 2 ■ 10~ 5 , D v = 10" 5 , (a) K = -0.4, r = 0.1 and (b) 
K = -0.8, r = 0.6. Control is switched on at t = 3000. 



spatio-temporal chaos has developed, control is switched on at t = 3000 and its effects on the time 
evolution are studied. Time duration of each run is 6000 in the regime of spatio-temporal chaos, 
and 7000 in the case of control of travelling waves. 

The results of numerical simulations of activator /inhibitor control are shown in Fig. |3^a) and 
(b). First of all, one can note a characteristic scaling of \K\ oc 1/t, which is quite natural bearing 
in mind that if the time delay is small, then the absolute value of the difference \u(t — r) — u(t)\ 
(or \v(t — t) — v(t)\) is small, and therefore a much higher value of the control strength \K\ is 
required to achieve a significant effect on the dynamics. In the case of activator control with K < 
or inhibitor control with K > 0, the only possible transition is from a spatio-temporally chaotic 
state to a non-trivial homogeneous state E\, as shown in Fig. |4j One should note that when the 
control strength is not sufficiently large, the system remains in the spatio-temporally chaotic state, 
as shown in Fig. Qa). The numerical boundary separating spatio-temporally chaotic state from a 
stable non-trivial homogeneous state shown in Figs. [4^a) andQb) coincides with an earlier found 
analytical expression for this boundary given by Eq. (|14[). 

In the opposite case, for activator control with K > or inhibitor control with K < 0, when 
the control strength \K\ and the time delay r are sufficiently small, the system remains in the state 
of spatio-temporal chaos (see Fig. (5^a)). For higher values of \K\t one observes stable mixed or 
Turing-Hopf modes shown in Fig. |5Jb) . These states are characterized by both spatial and temporal 
periodicity. Earlier work on co-dimension two Turing-Hopf states suggests that generically these 
are quite stable close to the Turing-Hopf boundary, and away from it they turn into either pure 
Hopf oscillations or Turing structures [26j [271 EH] • In the case of the Gray-Scott system this does 
not happen as the Turing instability boundary is completely enclosed inside the region of Hopf 
instability. 

For yet higher values of \K\t, isolated regions coalesce, and the system exhibits coarsening, as 
shown in Fig. p[c). Further increase in \K\t leads to a bistability between a homogeneous trivial 
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Figure 5: (Color online) Space-time plot in the case of inhibitor control, (a) Spatio-temporal chaos, 
(b) Mixed (Turing-Hopf) mode, (c) Coarsening, (d) Transition to travelling waves. Parameter 
values are a = 0.028, b = 0.053, D u = 2 • 10~ 5 , D v = 1(T 5 , (a) K = -0.05, r = 0.3, (b) K = -0.3, 
r = 0.35, (c) K = -0.4, r = 0.4 and (d) K = -0.6, r = 0.75. Control is switched on at t = 3000. 



state Eq or a localized pulse travelling with a constant velocity. This bistability here refers to a 
regime where for the same values of parameters and random initial conditions, the system will 
evolve either as shown in Fig. [5]^d), i.e. develop into a localized pulse, or it will develop in a manner 
qualitatively similar to the one shown in Fig. |4^b), except that the evolution leads to a uniform 
trivial steady state Eq. We also note that the time-delayed feedback of the form used in (10) is 
non-invasive only in the case of stabilization of a uniform steady state E\. In all other situations, 
although the feedback term does not vanish, its norm constitutes less than 10% compared to other 



terms in the system (10). 



Another possibility of control is the so-called diagonal control, which corresponds to the unity 
coupling matrix A. This type of control has been successfully used to stabilize unstable spatio- 
temporal breathing in a reaction-diffusion model for the double barrier resonant tunneling diode 
|15j . Numerical simulations suggest that for the Gray-Scott model this type of control does not 
qualitatively affect the dynamics as it leaves the system in a spatio-temporally chaotic state for 
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Figure 6: (Color online) Phase diagram for control of travelling waves, (a) Activator control, (b) 
Inhibitor control, (c) Phase diagram for non-diagonal control (15) of travelling waves. Colour 
denotes a specific final state: stable travelling wave (beige), stationary Turing pattern (cyan), 
bi-stability between travelling waves and a trivial steady state Eq (red). Parameter values are 
a = 0.022, b = 0.053, D u = 2 ■ 10" 5 , D v = 10~ 5 . 



any values of the control strength K and the time delay r. 

In the case of non-diagonal control that involves both species we introduce a variable phase 
into the coupling matrix as follows 



A 



cos <p 
— sin i 



sin ( 
cos ( 



(15) 



Such kind of control occurs, for example, in laser systems, where the optical phase can be used as 
an additional control parameter [291 E]- Figure [3^c) shows that for this type of control, the final 
state of the system is completely determined by the sign and size of non-diagonal control terms 
and shows no dependence on the values of the diagonal coupling terms. It is worth noting that 
in the case of negative feedback (—it < 4> < 0), i.e. when the roles of activator and inhibitor in 
the control term are reversed, one observes a range of possible dynamical regimes, such as stable 
Turing-Hopf modes, coarsening and bistability between travelling waves and a trivial steady state 
Eq. On the other hand, for positive feedback with < < 7r, the only possible transition is from a 
spatio-temporally chaotic state to a uniform non-trivial steady state E\. 

Beta and Mikhailov have previously shown for the complex Ginzburg-Landau equation that in 
the Benjamin-Feir regime, which gives rise to spatio-temporal chaos, the global TDFC is unable to 
stabilize uniform oscillations [9]. Our simulations suggest that a similar conclusion holds for local 
time-delayed feedback: while it is possible to recover homogeneous states, as well as stable mixed 
modes, it appears impossible to stabilize uniform oscillations in a reaction-diffusion system. 



4 Control of travelling waves 

In this section we consider a situation when without control the system supports stable travelling 
waves. By virtue of choosing periodic boundary conditions, if there is more than one pulse in the 
interval, all pulses are elastically reflected from each other and move without changes through the 
boundaries. 
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Figure 7: (Color online) Space-time plot of control of travelling waves, (a) Activator control: 
Development of a stationary Turing pattern, (b) Inhibitor control: Transition to travelling wave. 
Parameter values are a = 0.022, b = 0.053, D u = 2 ■ 10~ 5 , D v = 10~ 5 , (a) K = -0.5, r = 0.6 and 
(b) K = -0.8, t = 0.4. Control is switched on at t = 3000. 



As before, we start our analysis by considering single-species control of the form (11), and the 



results are shown in Fig. [6^a) and (b). In the case of activator control with K < or inhibitor 
control with K > 0, the stable travelling waves are eventually transformed into stable spatially 
periodic Turing patterns by means of wave-splitting, as shown in Fig. [7]^a). As has been noted by 
Petrov et al. [30j , this happens due to the fact that the medium becomes saturated with waves which 
are reflected from each other. Numerical simulations suggest that the higher the value of \K\t, 
the longer it will take for the wave splitting to result in a stationary spatially periodic pattern 
and that TDFC can stabilize unstable multi-spike steady states. This provides a mechanism of 
transformation from a stable propagating pulse to a stable Turing pattern. Finally, we note that 
the activator (K > 0) or inhibitor (K < 0) control shown in Fig. [6]^a) and (b), can only provide the 
bistability between a homogeneous trivial state Eq and travelling waves. For illustration purposes, 
we show in Figure [7tb) the transition to stable travelling pulses, however in the same parameter 
regime the system may also evolve into a homogeneous state Eq. 

Figure [6^c) shows the effects of mixed control on the dynamics of travelling waves with some 
similarity to the earlier picture of mixed control of spatio-temporal chaos (Fig. [3][c)). There is still 
almost no dependence on the strength and sign of activator control, but the region of stationary 
Turing patterns is noticeably larger than that of bistability between a homogeneous trivial steady 
state Eq and a travelling wave. As in the case of spatio-temporal chaos, the diagonal control of 
travelling waves does not qualitatively change the dynamics, and the system remains in the state 
of supporting elastic travelling waves. 

When compared to recent work on TDFC of travelling pulses in the FitzHugh-Nagumo model 
|17j . our results show both some similarities, and also certain differences. In both cases, the control 
domain of the inhibitor control is smaller than that of the activator control. It is worth noting 



that in [T7j non-diagonal control was not applied in the form of mixed control (15), involving both 



11 



(a) 



(b) 



0.8 



0.6 



0.4 



0.2 




0.8 



0.6 



0.4 



0.2 




-1 -0.5 0.5 1 



■1 -0.5 0.5 1 



Figure 8: (Color online) Phase diagram for control of travelling waves, (a) Activator itu-control. (b) 
Inhibitor uu-control. Colour denotes a specific final state: stable travelling wave (beige), stationary 
Turing pattern (cyan), bistability between travelling waves and a trivial steady state Eq (red). 
Parameter values are a = 0.022, b = 0.053, D u = 2 • 10 -5 , D v = 10~ 5 . 



activator and inhibitor variables, but was in the form of 
1 I J (uv — control), or A = 




±1 



(vu — control). 



Suppression of the traveling pulse in case of uu-control was only possible for the minus sign, while 
in case of -uu-control suppression was only found for the plus sign, following a strict symmetry with 
respect to the operation K — > —K. In some sense this is similar to Fig. |8]where the uniform steady 
state (red) can be stabilized only for negative uv feedback or positive vu feedback. 

The main difference, however, is that the non-diagonal control of travelling pulses in the Gray- 



Scott model in the form Eq. (15) is successful for both signs of the control strength K, even though 



for one sign of K it produces stationary Turing patterns and for another it leads to a homogeneous 
trivial steady state Eq. Finally, we note that tm-control is able to stabilize the uniform steady state 
or stationary Turing patterns for significantly smaller values of K and r as compared to ^-control. 
Also, Figure [8] shows that the roles of activator and inhibitor are reversed for the same values of K 
and r when compared to Fig. [6| 



5 Conclusions 

In this paper we have studied the effects of TDFC on the dynamics of the Gray-Scott model. When 
the original system is in the spatio-temporally chaotic state, an activator control with K < or 
inhibitor control with K > transforms it into a homogeneous non-trivial steady state E\. A 
much wider range of possible states can be reached for an activator control with K > or inhibitor 
control with K < 0, and these include the stable mixed Turing-Hopf modes, coarsening, as well as 
a regime of bistability between a stable pulse and a homogeneous trivial state Eq. Diagonal control 
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appears to be ineffective as it cannot control the spatio-temporal chaos. For mixed type of TDFC, 
the simulations suggest that the sign of non-diagonal terms in the coupling matrix plays a crucial 
role in determining the final state of the system. As in the case of global TDFC, the local control 
is unable to stabilize uniform periodic oscillations. 

When the TDFC is applied to travelling pulses of the Gray-Scott model, there are two possibil- 
ities for the final state. Activator control with K < or inhibitor control with K > transforms 
the system into a regime of wave splitting, and after the medium is saturated with the waves, it 
settles on a stationary spatially periodic Turing pattern. The higher is the value of \K\t, the longer 
it takes for the system to reach such a pattern. On the other hand, activator control with K > or 
inhibitor control with K < suppresses spatial inhomogeneities and provides bistability between 
a stable travelling pulse and a homogeneous trivial state, similar to the control of spatio-temporal 
chaos. 

The control schemes studied in this paper can be used for control of spatio-temporal chaos in 
other spatially extended systems. While stabilization of uniform periodic oscillations may remain 
elusive, local TDFC can provide a variety of interesting dynamical regimes, from stationary patterns 
to mixed modes and travelling localized pulses. Numerical simulations suggest that in many cases 
the control strength does not have to be very high, provided time-delay is large enough. This 
interplay between the control strength and time delay should prove beneficial for the development 
of experimental implementation of specific control schemes. 
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